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Cp ■ Abstract 

Curvature effects in Gamma-ray bursts (GRBs) have long been a source of consid- 
erable interest. In a collimated relativistic GRB jet, photons that are off-axis relative 
to the observer arrive at later times than on-axis photons and are also expected to 
be spectrally softer. In this work, we invoke a relatively simple kinematic two-shell 
collision model for a uniform jet profile and compare its predictions to GRB prompt- 
emission data for observations that have been attributed to curvature effects such as 
for the peak- flux-peak- frequency relation, i.e., the relation between the uFi, flux and 
the spectral peak, Epeak in the decay phase of a GRB pulse, and spectral lags. We 
present the case of the single-pulse Fermi GRB 110920, as a test for the predictions of 
the model against observations. 



1 Introduction 

Pulses in GRB light curves are thought to be produced by collisions between relativistic 
shells ejected from an active central engine (Rees and Meszaros 1994). The intercep- 
tion of a more slowly moving shell by a second shell that is ejected at a later time, 
but with faster speed, produces a shock that dissipates internal energy to accelerate 
the particles that emit the GRB radiation. This scenario is widely adopted in order 
to model pulses in GRB light curves (e.g., Daigne and Mochkovitch 1998; Zhang et al. 
2009). Studies of pulses are important to decide if GRB sources require engines that 
are long-lasting or impulsive, and to determine the likely radiation mechanism(s), with 
important implications for the nature of the central engine. 

Spectral lags, where low energy photons reach the observer at later times than high- 
energy photons, are seen in a significant fraction of GRBs. Cheng et al. (1995) were 
the first to analyze the spectral lag of GRBs, which they determined as the time delay 
between the peaks in the Burst and Transient Source Experiment (BATSE) Large Area 
Detector (LAD) channel 1 (25 - 50 keV) and channel 3 (100 - 300 keV) light curves. 
Since then, several authors have analyzed spectral lags in GRBs, while also extending 
these observations to the Swift and Fermi GRB samples (e.g., Norris et al. 1996; Norris, 
Marani & Bonneh 2000, Wu k Fenimore 2000; Chen et al. 2005; Ukwatta et al. 2010a) 
The leading model to explain the spectral lag is the curvature effect, i.e. the kinematic 
effect due to the observer looking at an increasingly off-axis annulus area relative to 
the line-of-sight (Salmonson 2000; loka and Nakamura 2001; Shen et al. 2005; Lu et al. 
2006). Softer low-energy radiation comes from the off-axis annulus area due to smaller 
Doppler factors. This radiation is also delayed at the observer end with respect to 
on-axis observation due to the geometric curvature of the shell. 

Both the existing models and observations suggest that a connection exists be- 
tween the observed hard-to-soft spectral evolution of GRB pulses and spectral lags. It 
is therefore important to understand the mechanism that produces this evolution. Ta- 
vani (1996) proposed that this hard-to-soft spectral evolution is caused by the variation 
of the average Lorentz factor of pre-accelerated particles and the strength of the local 
magnetic field at the GRB site as the synchrotron emission evolves within the burst. 
Liang (1997) proposed a physical model of hard-to-soft spectral evolution in which im- 
pulsively accelerated non-thermal leptons cool by saturated Compton up-scattering of 
soft photons. Kocevski &: Liang (2003) have pointed out that the nature of hard-to-soft 
spectral evolution is a difficult issue to address due to uncertainties involved with the 
microphysics of GRBs. They analyzed a sample of 19 GRBs, and found an empirical 
relation between the decay rate of the peak energy and the GRB lag. Ukwatta et al. 
(2010a) have analyzed a sample of 31 Swift GRBs with known red-shifts and deter- 
mined that the source-frame Ep^ak lies outside the energy band 100 - 250 keV for a 
majority of bursts, suggesting that spectral evolution may not be the dominant process 
causing the observed spectral lag. 

Borgonovo and Ryde (2001) studied the spectral evolution in the prompt-emission 
in GRBs by performing a time-resolved spectral analysis of BATSE single pulses and 
showed that in many cases the uFu flux at the Epeafc i Ie k^ f°^ each time segment 
was ex Epg^^ (hereafter referred to as the peak-flux-peak-frequency relation), with rj 



Table 1: Selected parameters for the generation of the model light curves unless otherwise stated. 
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ranging from ~ 0.6 to 3. The exponent r] was found to stay roughly constant for 
pulses within the same GRB. Dermer (2004) modeled and analyzed GRB pulses based 
on curvature effects with a Broken-Power-Law (BPL) intrinsic spectrum and showed 
that in the curvature limit, r] was equal to 3 for pulses with a wide range of temporal 
properties. He concluded on the basis of the Borgonovo and Ryde result that curvature 
effects do not play a prominent role in GRB prompt-emission. 

While several studies (Qin et al. 2003; 2004; Shen et al. 2005; Lu et al. 2005) have 
been performed to determine the role played by curvature effects, both in the spec- 
tral evolution of GRB prompt-emission as well as in producing spectral lags, many 
questions remain unanswered. An intriguing result connected to lag observations is 
the lag-luminosity anti-correlation (see Norris et al. 2000; Ukwatta et. al 2010a). This 
anti-correlation is further evidence for a connection between lags and the mechanism(s) 
that drive the GRB. 

In this paper, we investigate the effects of curvature by comparing the behavior 
of selective temporal and spectral characteristics of GRB prompt emission with the 
predictions of a simple two-shell collision model. The paper is organized as follows: 
the basic features of the model are presented in section 1.1, followed by a description of 
the sample selection criteria, analysis methodology, and a case study in section 2. The 
discussion of our main results is presented in section 3, followed by a brief summary of 
our conclusions in section 4. 

1.1 The Model 

We have used a particular representation of the internal shock model for our purposes 
(Dermer 2004). This model consists of a single two-shell collision event occurring 
at a radius ro from the source, generating a uniform spherical shell, with Lorentz 
factor r, that radiates for a co-moving time between tg and tg + At' with At' = 
ilA^ tvar/{^ + z). where t^ar is the observed variability time scale. The intrinsic- 
emission profile is assumed to be rectangular with instantaneous rise and decay phases. 
The co-moving width of this shell Ar' is assumed to remain constant during the period 
of illumination and given by Ar' = rj/^T c t^ar/i^ + z). The emission spectrum in the 
CO- moving frame may be described by any suitable spectral function such as, a Broken 
Power-Law (BPL), Band or Comptonized-Epgafe ; peaking at a co-moving photon energy 
E' ^ Q, = (1 + z) Epkfi/ 2 r where Ep^^o is the observer frame Ep^ at the start of the 
pulse. The spectral indices at E < Ep^^o and at E > Ep^^o in counts space are denoted 
a and /3, respectively. The curvature constraint requires that r < 2T'^c tvar/i^ + z) (see 
Dermer 2004 and references therein). The radius is thus written using the expression 
ro = 2r]r^'^c t^ar/i^ + z), with < r/j, < 1. The parameters rjt, ri/\ and rjr thus control 
the blast-wave duration, shell-thickness and radius of emission respectively. 

Unless otherwise stated, the selected parameters used for the numerical calculations 
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Figure 1: Normalized light curves obtained using se- 
lected parameters except: Thin-shell pulse, 77A — 0.1; 
and the Curvature pulse, »7a = ??t = 0.1. The case of 
rjt = rji:^ = rjr = 1 \s referred to as a Causal pulse (see 
Dermer 2004 for a detailed discussion of these three 
generic types of pulses). 
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Figure 2: Evolution of the spectral energy distribution 
due to curvature effects for the case of the causal pulse 
in Fig. 1 . In the declining phase of the pulse, the value 
of /e J. oc Epj., as shown by the red line. 
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Figure 3: Lag vs. Energy for selected parameters shell 
thickness t/a- Red: 77A = 1-0 cAt'; Green: t/a — 
0.5 cAt'; Blue: r]A = 0.1 cAt'. 
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Figure 4: Lag vs. Energy due to different intrin- 
sic spectral functions with selected parameters except 
Epfc.o = 200 keV. Here, Red: Broken-power-law with 
a ~ —2/3, /3 = —2.5, Green: Band function with 
a — —0.8, /3 = —2.25 and Blue: Comptonized-Epeafc 
with a = —0.8. 



of the light curves, spectra and spectral lags are shown in Table. 1. Fig. 1 shows the 
normalized, generic pulse shapes obtained using the selected parameters. Fig. 2 shows 
the evolution of the spectral energy distribution for the case of the curvature pulse 
shown in Fig. 1. The z^Fjy peak flux f £; j, oc E.^^ equality line is shown in the decay 
portion of the pulse. Dermer (2004) shows that this equality holds in the declining 
phase for all pulses. Fig. 3 shows the spectral lag as a function of energy for a case 
with selected parameters but with varying shell thickness, rj^. For a broken-power-law 
intrinsic spectrum, the lag shows a relatively well-defined "S-shaped" trend with no 
lags at energies below ~ O.SEp^ q aiid no lags for all energies above Ep^ q- Increasing or 
decreasing the value of Ep^^O) while keeping all other model parameters fixed, only shifts 
the entire curve in the direction of increase or decrease. Shen et al. (2005) studied 
the lags due to curvature effects using different intrinsic-emission profiles. They found 
that for an infinitesimal shell, a rectangular profile produces no lags. We find that the 
situation is quite different for the case of a finite shell thickness. Fig. 4 shows the lag 
as a function of energy for a case with selected parameters but with different intrinsic 
spectral functions (BPL, Band, and Comptonized-Epeafc)- It is seen that a change in 
the intrinsic spectrum has a significant effect in the evolution of the spectral lag with 
energy. In the case of a Band or Comptonized-Epeafc function, the lags are small for 
energies small compared to Epfc^O; and that the lags for a Comptonized-Epeafe function 
do not show a saturation energy. This is because the Comptonized-Epgafc function 
varies monotonically at all energies and does not have a well-defined break energy. 
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Figure 5: KRL pulse-fit for the light curve for GRB Figure 6: Light curve segments for 110920 with equal 
110920 with residuals. fluences. 



2 Sample Selection and Methodology. 

The model assumes that pulses in GRBs are produced by single two-shell-collision 
events. As a first step we therefore analyze either single-pulse GRBs, or GRBs with 
relatively simple light curves where the individual pulses within a multi-pulse structure 
in the light curve can be distinguished. In addition, we require that the GRB pulses 
be bright enough and of sufficient duration (the duration of the pulse is particularly 
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Figure 7: Comparable pulses generated using the BPL 
(Epeak = 300 keV, a = |, /3 = -i), Band (E^eafc = 
334 keV, a = -0.2, /3 = -2.65)and Comptonized- 
Epeafc (Epeafe = 280.1 kcV, a = —0.46) spectral func- 
tions. The light curves have been offset by 20 seconds 
for better viewing. 
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Figure 8: I'Fi' Flux vs. Epeak for the data from 
GRB 110920. The data were fit with the best- 
fit Comptonizcd-Epeafc function in the range 100-985 
keV. 
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Figure 9: I'Fi' Flux vs. Ep^ak from the model for the 
three intrinsic spectral functions described in the text 
with Blue: BPL; slope = 3.11 +/- 0.04, Pink: Band; 
slope = 3.39 +/- 0.10 and Black: Comptonized-Epeafc; 
slope = 2.57 +/- 0.01, for time segments identical to 
those used for the data. The flux scale has been offset 
for better viewing. 
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Figure 10: Lag vs. energy from the model for the 
pulses shown in Fig. 7 with Red: Data, Solid blue 
squares: Comptonized-Epeafc, Hollow black squares: 
Band function and Pink crosses: BPL. The Band 2 en- 
ergies are the mid-point of the energy in the second 
band. 



relevant to tests of the peak-flux - peak-frequency relation) so that we may obtain 
reliable results from our analyses. 

After identifying a potential candidate GRB, we pulse-fit the GRB light curves 
using a suitable pulse function (such as the Kocevski-Ryde-Liang (KRL) pulse func- 
tion: see Kocevski Ryde & Liang 2003 or the Norris pulse function: see Norris et al. 
2004) in multiple energy bands and then make a determination of the suitability of 
the GRB for our analyses. In order to confirm whether a given pulse (obtained from a 
suitable pulse-fit) is indeed a single pulse within statistics, we perform a wavelet based 
minimum-variability time-scale (MTS) extraction. The connection between pulses and 
the MTS has been established by MacLachlan et al. (2012). The best-fit pulse profile 
is then used as a representation of the light curve from the data. 

The time-integrated spectrum of the GRB is fit with a suitable function (Band, 
Comptonized-Epeafc etc.). The best-fit spectral function is used as the intrinsic emis- 
sion spectrum in the model. The model parameters are then varied to generate a light 
curve that best matches the best-fit pulse profile. The light curve is subdivided into 
time segments and an J^peak is extracted via a spectral fit for each time segment. The 
z/Fj, flux is extracted around ^peak- The model light curve is treated in an identical 
fashion as the data with regard to segments. Model fluxes and Epeafc's are extracted 
and the peak-fiux - peak- frequency relation is tested. Independent of this test, we 
also extract spectral lags in suitable, identical energy bands from the data and the 
model and compare the predicted and the observed spectral-lag-energy evolution. The 
spectral lags are extracted using the cross-correlation function (CCF) analysis method 
as described in Ukwatta et. al (2010). 



2.1 GRB 110920 - A Test Case 

The Fermi GRB 110920 is a single-pulse burst with a relatively long fast rise, expo- 
nential decay (FRED) structure with a Tqq of 170 =b 17 seconds. The best fit Band 
parameters (see McGlynn et al. 2012 for a detailed discussion on the properties of 
this GRB) for the time interval [Tq + 0.003, Tq + 52.737] (where Tq is the trigger time) 
were a = -0.20 ± 0.02, /3 = -2.65to:o9 and Epeak = 334 ± 5 keV with C-stat = 3206.5 
(485 d.o.f.). However, when a blackbody component was included in the fit, the C- 
stat was reduced to 2848.3 (483 d.o.f.). The peak energy of the Band component was 
shifted up to Ep^ak = 9'''8J]^2i ^^^ ^^^ the temperature of the blackbody was found 
to be kT = 61.3lQ'g keV. The low energy index a became shallower (—1.05 it 0.04). 
McGlynn et al. (2012) have attributed this blackbody component to the photospheric 
emission(see for instance: Ryde, 2005; Rees &: Meszaros, 2005). Since the redshift 
(z) is undetermined for this GRB, they have assumed a value of z = 2 and then used 
photospheric emission models to obtain a bulk Lorentz factor of ^^ 440. 

The persistence of the soft component at late times made initial attempts to pulse- fit 
the light curve rather difficult. A careful analysis, based partly on the results presented 
above, allowed us to separate the underlying pulse structure from the overall structure 
of the light curve. Fig. 5 shows the full light curve for the GRB in the energy range 100 



- 985 keV along with the best-fit KRL pulse function. The energy range was so chosen 
to avoid contamination from the soft component (below 100 keV). Once this was done, 
the parameters obtained from the KRL pulse-fits in different energy bands showed a 
consistent pattern in terms of their values, as well as the degree of correlation between 
them. The MTS we extracted for this light curve is about 1.3 seconds, thus suggesting 
that this GRB is most likely a single-pulse GRB with no additional structures. 

The best-fit KRL function was then used as the representation of the light curve 
to be matched by the model. We assumed a Bulk Lorentz factor T of 440 and z = 
2 (identical to McGlynn et al. 2012). With these parameters, we extract a radius, rg 
of about 10^^ cm from the model. Fig. 6 shows the segmentation of the light curve 
into equal-fluence segments. The time-resolved energy-spectrum of the pulse was fit 
with a Band function (identical to McGlynn et al 2012, in the energy range, 8 - 985 
keV) and a Comptonized-Epeafc function (in the range 100-985 keV) and these, along 
with a theoretical BPL function (originally used by Dermer 2004) were used as in- 
trinsic spectra for the model. With a suitable choice of shell-collision parameters, we 
extracted best-fit pulses using these three spectral functions (see Fig. 7 for the pulses 
as well as details of the corresponding spectral parameters). Figs. 8 and 9 show our 
results for the peak-flux-peak-frequency for the data with the best-fit Comptonized- 
Epeafc function and the model using the three different intrinsic spectra. Fig. 10 shows 
the spectral lags for the light curves of Fig. 7. 



3 Discussion 

The above results show that an internal shock model with an intrinsic spectrum iden- 
tical to that used to fit the data (the Comptonized-Epeafc and Band functions), repro- 
duces the observed pulse profile as well as the observed spectral lags. It was difficult to 
determine if there was a saturation energy present in the lag-energy-evolution (Fig. 10) 
as there were insufficient counts to extract lags in higher energy bands. As noted earlier, 
a finite shell thickness can account for observed lags even for the case of a rectangular 
intrinsic emission profile. A second Band function component with a peak shifted to 1 
MeV (when a blackbody component is included in the fit) would imply a 1 MeV break 
energy in the lag-energy plot, and would also predict no lags (or small lags in the case 
of the Comptonized-Epeafc) below ~300 keV. This does not match the observations. It 
can be seen from Fig. 9 that the exponent in the peak-flux-peak-frequency relation is 
close to 3 in all cases for the model pulses shown in Fig. 7. This conflrms the observa- 
tions of Dermer (2004) that the exponent at times after the peak of the light curve is 
close to 3 even with different choices of intrinsic spectra. However, this does not match 
the exponent obtained from the data (Fig. 8) . This further supports the conclusions of 
Dermer (2004) that curvature alone cannot account for the spectral evolution seen in 
GRB prompt-emission. While a connection may exist between the observed spectral 
evolution and the spectral lags, it appears that curvature effects alone cannot describe 
this connection and additional emission mechanisms may be needed. We note in pass- 
ing that Guirec et al. (2012) have shown that the inclusion of a blackbody component 
in the spectral flts leads to an improved description of the spectral properties of a 
number of GRBs. 



4 Conclusions 

We have used a simple two-shell collision model to investigate curvature effects in the 
prompt emission of GRBs. In addition to the BPL for the intrinsic-emission spec- 
trum used in Dermer (2004), we have also examined the effects of using other emission 
spectra such as the Band, and the Comptonized-Epeafc functions. In this study, we fo- 
cused primarily on the peak-flux - peak-frequency relation which is predicted to have 
a power-law behavior with an exponent of 3, and the energy evolution of spectral lags. 
We compare our model results with the results of similar models in the literature and 
also present a test case study of GRB 110920. 

We summarize our main findings as follows: 

• While we agree with Shen et al (2005) that an infinitesimal shell produces no 
discernible spectral lag using a rectangular emission-pulse-profile, we find the 
situation to be different for a shell of finite thickness i.e., a finite spectral lag can 
be produced even with a rectangular pulse profile; 

• The spectral lag evolution as a function of energy is quite sensitive to the type of 
intrinsic spectrum. For example, the Comptonized-Epeafc rnodel does not appear 
to exhibit a saturation energy at which the spectral lags reach a plateau phase as 
in the case of the Band and the broken-power-law functions. We agree with Shen 
at al (2005) that the spectral lags seem to approach a maximum when Epeafc is 
near the high-energy channel used in extracting the lag. Most likely this simply 
reflects the break energy present in the assumed intrinsic spectrum (i.e.. Band 
and BPL); 

• All intrinsic spectral models tested exhibit the peak-flux - peak-frequency rela- 
tion although with exponents that differ from the predicted exponent of 3. The 
significance of this discrepancy is not clear at this stage and warrants further 
investigation; 

• The peak-flux - peak-frequency test for GRB 110920 yields an exponent of 1.64 
+/- 0.012 compared to the theoretical one of 2.57 +/- 0.01 (with Comptonized- 
Epeafc as the intrinsic spectrum). We consider this discrepancy to be significant 
and the result to be in disagreement with the prediction based purely on effects 
of curvature. Similar conclusions were reached by Dermer (2004) and Borgonovo 
and Ryde (2001); 

• Relatively good agreement is obtained with all intrinsic spectral models for the 
spectral lag versus energy for GRB 110920. This is somewhat surprising given 
the result stated in the previous bullet. It is possible that this result is due to 
the interplay between model parameters such as shell thickness, variability time 
scale, the energy evolution of Epeafc and the Lorentz factor. The investigation of 
the dependencies of these various parameters is ongoing. 

We note that the role of the reported soft component of the light curve for GRB 
110920 has not been fully investigated in this study and is worth pursuing, particularly 
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with regard to the behavior of the peak-flux-peak-frequency relation. Finally, we note 
that these studies are being extended to a larger sample of GRBs. 
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